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ABSTRACT 

Meson-exchange current (MEC) contributions to the parity- violating (PV) asymmetry 
for elastic scattering of polarized electrons from 4 He are calculated over a range of momen- 
tum transfer using Monte Carlo methods and a variational 4 He ground state wavefunction. 
The results indicate that MEC's generate a negligible contribution to the asymmetry at 
low-|cf|, where a determination of the nucleon's mean square strangeness radius could be 
carried out at CEBAF. At larger values of momentum transfer - beyond the first diffrac- 
tion minimum - two-body corrections from the p-n "strangeness charge" operator enter the 
asymmetry at a potentially observable level, even in the limit of vanishing strange-quark 
matrix elements of the nucleon. For purposes of constraining the nucleon's strangeness 
electric form factor, theoretical uncertainties associated with these MEC contributions do 
not appear to impose serious limitations. 
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I. Introduction. One objective of the CEBAF physics program is to probe the strange- 
quark "content" of the nucleon with parity-violating (PV) electron scattering. As dis- 
cussed elsewhere in the literature [1-7], PV electron scattering at low-to-intermediate en- 
ergies is particularly suited to the study of strange-quark vector current matrix elements, 
(H\s'Jns\H), where if is a hadron. In the case where the target is a nucleon (\H) = \p) or 
\n)), this matrix element can parameterized by two form factors, Ge(Q 2 ) and Gm\q 2 ), 
the strangeness electric and magnetic form factors, respectively. Extractions of (N\ss\N), 
the nucleon's strange-quark scalar density, from tt—N scattering [8,9], as well as determina- 
tions of the strange-quark axial vector matrix element, (N | ^7^75 s\N), from elastic VnP/v^p 
scattering [10-12] and measurements of the gi sum [13-15], suggest that the strange-quark 
"sea" plays a more important role in the low-energy properties of the nucleon than one 
might expect based on the success of valence quark models. Measurements of (iV|s7 M s|iV) 
would provide an additional window on the sea-quark structure of the nucleon. Model 

(s) (s) 9 1 

estimates of G E and G M at low-|Q | span a wide spectrum in both magnitude and sign 
[16-21]. It is therefore of interest to extract the strangeness form factors at a level needed 
to distinguish among model calculations and their attendant physical pictures. 

To this end, use of a proton target would not be sufficient. The presence of several 
poorly-constrained form factors in the PV elastic 1 H(e, e) asymmetry, as well as theoretical 
uncertainties associated with axial vector radiative corrections, limit the precision with 
which G^ and G$ could be determined from the proton alone [1,2]. The use of A > 1 
targets in conjuction with the proton offers the possibility of imposing more stringent limits 
on the nucleon's s-quark vector current matrix elements [1,2,22] than could be obtained 
with a proton target only. In this regard, the (J n ,T) = (0 + ,0) nuclei, such as 4 He, 
constitute an attractive case, since the ground states of such nuclei can support matrix 
elements of only one operator - the isoscalar Coulomb operator [1,2,22,23]. In the one- 
body approximation to this operator, the nuclear wavefunction dependence of the Coulomb 
matrix elements effectively cancels out from the PV asymmetry for such nuclei, leaving 
only a sensitivity to Standard Model couplings and single nucleon form factors (e.g., G^). 
Two approved CEBAF experiments rely on this feature of A LR (0 + ,0), the PV left-right 
asymmetry [24,25]. The proper interpretation of A LR (0 + ,0) requires that one understand 
the importance of many-body corrections to the one-body asymmetry. Meson-exchange 
currents (MEC's) constitute one class of such many-body effects. In previous work [26], 
we computed MEC contributions to the 4 He mean-square "strangeness radius", which 
generates the leading s-quark contribution to A LR ( 4 He) at low-|cf|. The results of that 
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calculation, peformed with a simple 4 He shell model wavefunction and phenomenological 
two-body correlation function, indicate that the 4 He strangeness radius is dominated by 
strange-quarks inside the nucleon. 

In the present work, we extend the calculation of Ref. [26] using a 4 He variational 
wavefunction obtained from realistic interactions and computing the asymmetry over the 
full range of momentum transfer germane to the future CEBAF experiments. Our results 
indicate that the 4 He strangeness radius is two orders of magnitude more sensitive to the 
nucleon's strangeness radius than to two-body contributions. At the higher \q\ of experi- 
ment [24], the situation is more complex. Even if the nucleon matrix element (A^|s7^s|A^) 
were to vanish, the PV asymmetry would still receive a non-negligible contribution from 
non-nucleonic s-quark matrix elements. In particular, the p — tv strangeness transition 
charge operator generates nearly a 15% contribution to the asymmetry at the kinematics 
of the experiment [24] . In this case, an experiment like that of Ref. [24] would be significant 
in two respects. First, it would be interesting to measure a non-negligible strange-quark 
matrix element in a strongly-interacting, non-strange system, regardless of the dynamical 
origin of that matrix element. Second, the only other observable with significant sensitivity 
to the p — 7r MEC is the B form factor of the deuteron [27]. If, however, Ge an d Gm' are 
non-zero, the level of theoretical uncertainty associated with the present MEC calculation 

(s) 

does not appear to be large enough to significantly weaken the possible constraints on G E 
which a measurement of A LR ( 4 He) could provide. 

In the remainder of the paper we provide details of the calculations leading to these 
conclusions. Section II gives our formalism, including expressions for the operators used. 
In section III, we treat the computation of the 4 He matrix elements of these operators, 
considering first the simple case of a shell model ground state and subsequently turning to 
the Variational Monte Carlo (VMC) approach. In section IV we discuss our results, includ- 
ing implications for the interpretation of A LR ( 4 He) and studies of nucleonic strangeness. 
Technical details may be found in the Appendix. 

II. Formalism. The PV left-right asymmetry for scattering of polarized electrons from 
a nuclear target depends on the interference of the electromagnetic (EM) and PV weak 
neutral current (NC) amplitudes, M EM and M££, as 

^ 2Re Ml M MZ 

\M EM \- ' (1) 
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where \M EM \ » \M™\ at low energies. The amplitude M PV C is proportional to the sum 
of two terms, 

fH JNC , pfl5 JNC (r>\ 

where is the electron's vector (axial vector) neutral current and ^(5) is the nucleon 
or nuclear matrix element of the hadronic vector (axial vector) NC. One may rewrite A LR 
in terms of quantities which set the scale of the asymmetry and a ratio of nuclear response 
functions [1,2] 

where G M is the Fermi constant measured in muon-decay, a is the EM fine structure 
constant, and Q 2 = u 2 — q 2 with u> and q = \q\ being the energy and magnitude of three- 
momentum transfer to the target. The response functions appearing in the ratio of Eq. (3) 
may be written as 

F 2 =v L R L +v T R T (4a) 
W PV =v L WZ v + v T WT v + v T ,W? A , (4b) 

where v L , v T and v T i are leptonic kinematic factors; R L and R T are the usual longitudinal 
and transverse EM response functions; and and Wy' A are analogous PV response 

functions involving products of the hadronic EM and vector NC ( u AV n ) or axial vector 
NC ("VA") [1,2]. 

In this work, we follow the approach taken in Refs. [1-7] and keep only the three 
lightest quarks in the hadronic current. In this case, one has for the two vector currents 

jem = j^( T = 1) + J° M {T = 0) (5a) 
J»° = C =1 J^ C ( T = 1) + ^1 =G J; M {T = 0) + tfh^s (5b) 

where the J^ M (T) are the isovector (T = 1) and isoscalar (T = 0) EM currents and the 
£v are couplings determined by the Standard Model [1,2,22]. A decomposition of J^ 5 C 
analogous to that of Eq. (5b) but involving the SU(3) octet of axial currents and s^^ 75s 
may also be made [1,2,7]. Since the 4 He ground state supports no axial vector matrix 
element, however, we do not consider J^ 5 C further in this work. 

In the limit that the 4 He ground state is an eigenstate of isospin, the "hadronic ratio" 
for this target is 



F 2 2 ■•" s " F^°(q) 
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Here, V3^ =0 = -4sin 2 6 W and ^ 0) = -1 at tree level in the Standard Model [1,2,22]. 
The form factors are given by 

F^(q) = (0+0\M$(q)\0+0) (7a) 
M${q) = j d 3 xjo(qx)Y o(n x )p^(x) (7b) 

= ±-f dn q Y 00 (n q )p^(q) , (7c) 

where x = \x\ and p^ a \x) (p^ a \q)) denotes the co-ordinate-space (momentum-space) 
charge (p = 0) component of either the isoscalar EM current ((a) — > T = 0) or strange 
quark current ((a) — > (s)). Matrix elements of the Coulomb operator are simply related 
to the elastic charge form factor as 

4%) = (0+0|p(<f)|0 + 0> = 2V5FF^(g) . (7d) 

One observes from Eq. (6) that were the nuclear matrix elements of Mqq (q) to vanish, the 
asymmetry would be nominally independent of the details of the nuclear wavefunction.f 
The reason is that (i) in the absence of strangeness, the hadronic isoscalar EM and isoscalar 
NC currents are identical, up to the overall electroweak coupling, \/3£v =( \ (ii) isovector 
matrix elements vanish if the 4 He ground state is assumed to be a pure T = state, and 
(iii) a spin-0 ground state cannot support axial vector matrix elements. 

One-body operators 

Expressions for the one-body charge operators may be obtained starting from Lorentz- 
covariant forms of the single- nucleon vector current matrix element: 



(N(p')\V,(0)\N(p)) = U(p') 



2m N 



U(p) , (8) 



where Fi and F 2 are the standard Dirac and Pauli form factors of the nucleon, U(p) and 
U(jp') are nucleon spinors corresponding to nucleon states \N(p)) and \N(p')), respectively, 
and Vfj,(x) is any one of the vector currents of interest (isoscalar EM or strangeness). 



f Apart from contributions from nuclear dispersion corrections; see, e.g. Refs. [1,2,22]. 
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Expanding the right side of Eq. (8) in powers of p/m N , transforming to co-ordinate space, 
and summing over all nucleons gives for the \i = component 



fc=l 
)2 IA^2 



8m 2 , 



G { E a) (T)-2G^\r)}a k -qxP k 



(9) 



where r = —Q 2 /Am 2 N = g 2 /4m 2 for elastic scattering in the Breit frame, P k = p k + P k i 
and 



G { E a) 

Ma) 

M 



F (a) _ rF {a) 
F (a) + p(a) 



(10a) 

(10b) 



are the Sachs electric (10a) and magnetic (10b) form factors [28]. In arriving at the 
expression in Eq. (9), we have used the spinor normalization of Ref. [29]. Had we followed 
the convention of Ref. [30] , the charge operator would have contained an additional term 
A k inside the square brackets given by 

1 



A fc = 



Am 2 ., 



Pk + P'Vj (1 + r ) _1 [G { E a \r)+rG[f(r) 



(11) 



Following the convention in Refs. [1-4], we parameterize the momentum-dependence 
of the one-body form factors as 



GI(t) = G°(t) 
<? m {t) = » p G d v {t) 



(12a) 
(12b) 



GI{t) = -^tG d v {t)Ut) 



(12c) 
(12d) 



with 



G { E s \r) 
G ( °\r) 



G° v {r) 

is 



Ps tG d v (t)Ut) 
HsG d v {t) 



(l + Xnr)- 1 
(l + A^V)- 1 



(12e) 
(12f) 



(13a) 
(13b) 
(13c) 
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and 



GlT° = \ [<? BtM + Gl M ] (14a) 
GT M = \ [Gl M - Gl M ] . (14b) 

Numerically, one has \i v = 2.79, [i n = —1.91, X° = 4.97, and A n = 5.6. The rationale for 
adopting this parameterization is discussed more fully in Refs. [1,2]. The parameters fi s and 
p s , which define the strangeness magnetic moment and strangeness radius, respectively, as 
well as which governs the next-to-leading Q 2 behavior of G^\ are presently unknown. 
One goal of the SAMPLE experiment [6] and up-coming CEBAF experiments [24,25,31] 
is to place limits on these parameters. 

The one-body contribution to the Coulomb multipole operator, obtained by substi- 
tuting the expression for the charge operator of Eq. (9) into Eq. (7b), is 

^o\^) [1] = ^^1^=^) (15) 



+ 



G^\t)-2G^\t) 



J 2m N m N Xk 

where we have assumed qo = so that q 2 = Am 2 N r and where Lk is the orbital angular 
momentum of the /c-th nucleon. Note that in the limit that the 4 He ground state consists 
of nucleons in S-states only, the spin-orbit operator in M^ 1 (q)^ will not contribute to 
Fqq (q) . In this case, the Coulomb matrix elements for the isoscalar EM and strangeness 
charge operators are identical, apart from the single nucleon form factors, rendering their 
ratio independent of nuclear structure: 



G { E s \r) 



S— waves E 



GT°{r) 



(16) 



However, the presence of a significant D-wave component (the associated probability is 
about 16% for the variational He wavefunction discussed below) implies some level of 
structure-dependence in the one-body form factor ratio of Eq. (16). For most values of 
momentum transfer, the magnitude of this structure-dependence is negligible (see the spin- 
orbit contributions in Figs. 2b and 3b). 



Two-body operators 

In the one boson-exchange (OBE) approximation, the leading two-body MEC correc- 
tions to the one-body result of Eq. (16) are generated by the processes in Fig. 1. The 
7r-exchange and vector meson-exchange "pair currents" (Figs. la,b) are familiar from pre- 
vious work on MEC's [32-36], as is the pseudoscalar- vector meson "transition current" of 
Fig. lc. In each case, the isoscalar EM and strangeness two-body currents have the same 
structure, apart from the form factors appearing at the NN creation/annihilation vertex 
and V — n transition vertex. The 4 He elastic from factors receive no contribution from 
processes in which a virtual 7 or Z° couples to an exchanged pseudoscalar or vector me- 
son. The reason is that matrix elements of the form (M\V^(0)\M') must vanish in order 
to respect G-parity invariance when \M') and \M) are identical meson states (apart from 
momenta) and when is either J^ M (T = 0) or s7 M s. Moreover, one has no contribution 
from an u> — it transition current since currents which are strong isoscalar operators can- 
not induce such an isospin-changing transition. We have not included contributions from 
isobar currents, since the lightest nucleon resonance accessible with an isoscalar current is 
the AT(1440). We assume contributions from the associated current are suppressed by the 
large mass difference between this state and the nucleon. 

We derive two-body charge operators by computing the covariant momentum-space 
Feynman amplitudes associated with the diagrams in Fig. 1, peforming the standard non- 
relativistic reduction, and transforming to co-ordinate space. We take the meson-nucleon 
couplings from the conventional low-energy effective Lagrangians: 



jCnntt = -^ LI1IL 'ipN(x)p7f(x) ■ tiI> n (x) (17a) 
2m N 

£nn p = 9 p nn^n{x) 

and 

= d» + ieQ EM A» + igQ w Z fl , (17d) 

where V'jV is a nucleon field, 7r a , p®, and cu M are the pion, rho- meson, and omega- meson 
fields, respectively, "a" is an isospin index, g is the semi-weak coupling, Q EM and Q w are 



7 M + 



Kp 

2m, 



Pfi ■ T1p N {x) 



(17b) 



2m, 



u^ip N {x) 



(17c) 
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the EM and weak NC charge operators, and and are the photon and Z° fields, 
respectively. We take the couplings appearing in Eq. (17) to have the values g TiVJV = 13.6, 
9 p nn = 2.6, g UNN = 14.6, k p = 6.6, and k u = —0.12 [37]. Momentum-space matrix 
elements of the operators in Eq. (17) have the same structure as the effective Lagrangians, 
but with the nucleon fields replaced by plane wave spinors, the derivatives replaced by 
where /c M is the momentum of the outgoing meson, and the vector boson fields replaced by 
the corresponding polarization vectors, £ M . For the p — 7T transition current matrix element 
one has 

(7r%)|^°>(0)|p c (fci,e)) = - 9 ^ {Q2) 5 bc e^ aP k^ , (18) 

Vtlp 

where as usual "a" denotes either the EM or strange-quark current [38]. In the case of 
the former, the value of the transition form factor at the photon point is known to be 
gp=°(Q 2 = 0) = <7 p7r7 = 0.56 [39], while the Q 2 -dependence may be modelled using w-pole 
dominance: 

9 T p ;°(Q 2 )=9 P «A 1 -Q 2 / m ir 1 ■ (19) 

In the case where = «7 M s, one may follow a similar approach and assume <^>-meson 
dominance, which is reasonable since the 4> is almost pure ss: 

9 { ;KQ 2 ) = 9 P ,s{l-Q 2 /mlY 1 • (20) 

The measured rates for — > pir and <j> — * £ + £~ (£ is a charged lepton) can be used to 
estimate the value of this form factor at Q 2 = to be \g p7rs \ = 0.26 [40]. 

Before proceeding, we touch on one issue associated with the vector meson pair cur- 
rent operators. These operators are derived by keeping only the negative-energy pole of 
the nucleon propagator, as shown in Fig. 1. The resulting two-body nuclear matrix ele- 
ment is thus distinct from the matrix element containing the positive energy pole, which 
contributes via the full nuclear Green's function in time-ordered perturbation theory: 



E 



(f\J P \n){n\H NUC \i) (f\H NUC \n){n\J^\i) 
Ei~ E n + is Ei- E n + is 



(21) 



Here, H NUC is the full nuclear Hamiltonian and (i, f, n) denote initial, final, and inter- 
mediate nuclear states, respectively. Following this prescription leads one to a two-body 
pair-current operator having the same form as given in Ref. [36]. It has been argued, 
however, that one must include an additional retardation contribution arising from the 
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positive-energy pole in the nucleon propagator whose residue contains a dependence on 
the energy transfer between the two nucleons. Inclusion of this additional term results in 
the form for the pair-current charge operator given in Refs. [32,33]. Rather than attempt- 
ing to choose between these two approaches, we compute Fqq in two ways - once using 
each of these two prescriptions - in order to determine the impact of this choice. As we 
note in Section IV, the vector-meson exchange contributions to the 4 He form factors are 
sufficiently small in comparison with other contributions that the impact of this choice in 
the value for Fqq is insignificant. 

The momentum-space charge operators for the pair currents are 



p(pi,Pi,P2,P2> £)pi ] onic = (2nf5(k 1 + k 2 -q) 

1 



9-KNN 

4ml 



k% + ml 



<7i ■ qa 2 • k 2 + (1 <-> 2) 



(22a) 



where ki = p{ — fy, i = 1, 2, and 



p{puPi,P2^ q) [ v ] J;L = (2nf5(k 1 + k 2 - q) 



9vNN 

4ml 



x 



1 



k\ +ml 
-iai xq - (p 2 +p 2 ) 



(1 + k v ) (q- k 2 + a i x q- a 2 x k 2 j 
+ (1 <->2) 



(22b) 



excluding the retardation correction or 



m,P(,P2,P 2 ;q) [ yl { X = (2nf5(h + k 2 -q) 



(1 + k v ) 2 F?\t) + (1 + k v )F?>(t) ) S 1 x q- a 2 x k 2 



k\ + m^, 



9vNN 

4ml 



Tv(l,2) 



+ Gm\t) ( (1 + K v )q ■ k 2 - iai x q ■ (p 2 + p 2 ) 



+ (!<-> 2) 



(22c) 



including the retardation term, with 



f v (l,2) 



T\ ■ T 2 , V = p 

1 , V = iv 



(23) 



10 



and where "(i)" indicates either the isoscalar EM or strange-quark charge operators. We 
have not included the isovector parts of the charge operators. For the transition operators, 
one has 



p(pi,pl,P2,P 2 ; Q) [ ol = ^(27^) 3 ^(fc 1 + k 2 -q) 



Am p m 2 N 



(24) 



x 



-(Ti • ki 



(kj + ml)(kl + ml) 
- i(l + K p )(ki x k 2 ) ■ (k 2 x a 2 ) 



(pi + p[) ■ (h x k 2 ) 
+ (!<-> 2) 



Expressions for the co-ordinate space forms of the two-body charge operators, 
p(xi, x[, X2, x 2 ] (?) , as well as for their Coulomb multipole projections, Moo(<?) > are 
somewhat involved and may be found in the Appendix. For purposes of discussion, it 
is useful to consider the leading-g behavior of the two-body Coulomb operators (shown 
in Eqs. (A. 10) of the Appendix), since their matrix elements contribute to the 4 He EM 
and strangeness radii. From the \ow-q expressions for the two-body Coulomb operators, 
we observe that they vanish at least as rapidly as q 2 for small q. The operators must 
vanish at q 2 = 0, since the two-body operators cannot change the overall charge (EM or 
strangeness) of the 4 He nucleus. In the case of strangeness, the entire nuclear form factor 
F$ must vanish at q 2 = 0, since the nucleus has no net strangeness. Thus, in analogy 
with the single nucleon case, we define a nuclear strangeness radius as 



p s [nuc] = 2a/7t 



dF. 



(s) 



CO 



dr 



(25) 



T = 



Under this definition, p s [nuc] = Ap s in the one-body limit neglecting the spin-orbit contri- 
bution. From the expressions in Eq. (A. 10), we note that the pionic operator (Eq. (A. 10a)) 



contributes to Fqq at C(i? 4 ), since this operator is proportional to q 2 F\ a> and since F[ 
vanishes as q 2 for small q. Consequently, the longest-range MEC does not contribute to 
the nuclear strangeness radius. For the same reason, the retardation correction to the 
vector meson pair current operator (Eq. (22c) and Ref. [33]) also does not contribute to 
p s [nuc], since this correction is proportional to tF^ s \t). As a result, the low-g behavior of 
the vector meson contribution to F^q is independent of the choice of approach discussed 
above. This choice takes on relevance only at larger values of momentum-transfer, where 
the terms proportional to tF^ s \t) are non-negligible. 



(a) 
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III. Calculation of 4 He Matrix Elements. Although the object of this paper is 
to report on a calculation of F^° and using state of the art wavefunctions, we 

first summarize a simpler calculation of the 4 He strangeness radius using a shell model 
ground state with harmonic oscillator wavefunctions. This simpler treatment allows for 
an analytical computation and serves to guide one's intuition when interpreting results 
obtained with more sophisticated methods. The results of the shell model calculation were 
reported previously [26], and we provide more details in the Appendix of the present paper. 

Shell Model Calculation 

In the simplest shell model description of 4 He, the ground state consists of a single 
configuration: four nucleons in the lsi/ 2 state. Numerical results using more realistic 
wavefunctions, such as the variational wavefunction described below, suggest that the level 
of configuration mixing is at least 15% . Within the S-state approximation, we compute 
the leading-g behavior of Fqq using harmonic oscillator single-particle wavefunctions with 
an oscillator parameter b = 1.2 fm, obtained from fits to the data on Fqq° [2]. Analytic 
expressions for the nuclear matrix elements appear in the Appendix, and our results give 

F«(r - 0) = ^rp/He] (26) 

= T (XlPs + A 2a ^s + A 2 bfi f p7rs) , 

where the terms containing Ai and A2 a ,b give the one- and two-body contributions, respec- 
tively. The one-body term is nuclear structure-independent, since the leading (/-dependence 

(s) 

of the one-body strangeness Coulomb operator is given by G E times an operator which 
counts the number of nucleons (see Eq. (15)). The two-body term \2&p s arises from 
the vector meson pair currents, while the term \2h9pirs is generated by the p-n transition 
current. Numerically, in the limit of point meson-nucleon vertices (A M — > 00), we obtain 
Ai rs 1.13, A 2a ~ —0.05, and A 2 b ~ —0.02 after including a phenomenological NN anti- 
correlation function in the two-body matrix elements. We expect that the values of the 
A2a,b fc> r finite A M should be smaller in magnitude than those quoted, which we take to 
give an upper bound on the scale of two-body contributions. These results imply, then, 
that p s [ 4 He] is at least a factor of 20 more sensitive to the nucleon's strangeness radius 
than to two-body strangeness currents. 

We note in passing that had we not accounted for short range NN repulsion, the 
vector meson contribution would have been a factor two larger in magnitude and the p-n 
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term matrix element would have been a factor of ten larger. The reason for the large 
suppression of the p-ir term due to short-range repulsion can be seen from the structure of 
the momentum-space p-7r charge operator in Eq. (24). At leading-order in q, the Coulomb 
projection of this operator has the form 



q 2 



(k\ + ml)(kl + m 2 ) 



+ 0(q 4 ) 



(27) 



{m 2 - ml) 



O 



k\ + ml k 2 + 



rat 



+ 0(q 4 ) 



where O is an operator dependent on and k\. Nuclear matrix elements of the full 
operator in Eq. (27) thus depend on the difference of matrix elements of two operators, 
A(m n ) and A(m p ), whose ranges are set by ra^ and m p , respectively. In the absence 
of short-range anti-correlations, one has 2(g.s. ^(ra^)!! g.s.) ~ (g.s. ||A(ra p )|| g.s.). The 
impact of short range repulsion is to reduce the p-meson term (g.s. || A(m p ) \\ g.s.) by about 
a factor of two, while leaving the matrix element of the pionic operator, whose range is 
much larger than the radius of the repulsive core, relatively unchanged. Consequently, the 
degree of cancellation between the two pieces is greatly enhanced, leading to the factor of 
ten reduction in A2t>, as compared with the less significant impact on the magnitude of the 
purely vector meson matrix elements, A2a- 

Variational Monte Carlo Calculation 

The 4 He variational wavefunction used in the present work is obtained by minimizing 
a realistic Hamiltonian with the Argonne V14 two-nucleon [41] and Urbana-VIII three- 
nucleon [42] interaction models. It has the symmetrized product form given by [42]: 



i*>=[i+ y: u5n[siki+ua)]\*3> 

i<j<k i<j 

Here S is the symmetrizer, and > is a Jastrow wavefunction 

i<j 



(28) 



(29) 



where A is the antisymmetrizer acting on the spin-isospin states of the four nucleons. The 
two-body correlation operator Uij is taken to be 



P=T,a,(TT,t,tT 



(30) 
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with 

Ofj = Ti • Tj , (Jj • (Tj , (Ti • (TjTi • Tj , Sij , SV/Ti • f j ; p = T, <T, (7T, t, £t . (31) 

The three-body correlation operator U™ 1 is simply related to the three-nucleon interac- 
tion present in the Hamiltonian, and has a correspondingly complex operator dependence. 
The correlation functions f c (r) and u p (r) as well as the additional parameters present in 
U™ 1 are determined variationally with the methods discussed in detail in ref. [42] . 

The 4 He binding energy and charge radius calculated with the above wavefunction 
have errors of ~ 4% when compared to exact Green's Function Monte Carlo (GFMC) 
results for the same Hamiltonian [42,43] (we note that the GFMC results reproduce the 
empirical values). This wavefunction also produces a charge form factor that is in good 
agreement with the exact GFMC predictions and the experimental data over a wide range 
of momentum transfers [42]. Because of the relatively strong tensor component in the 
Argonne V14 the D-state probability has the rather large value of 16%. 

The charge and strangeness form factors are given by the expectation values 

4 a) = 2V^4o(<7) =< y;q\p (a \m > , (32) 

where \^;q > denotes the ground state wavefunction recoiling with momentum g, and 
p^ a \q) are the r-space representations of the charge and strangeness operators listed in 
the appendix. The above expectation value is computed, without any approximation, by 
Monte Carlo integration. The wavefunction is written as a vector in the spin-isospin space 
of the A-nucleons for any given spatial configuration R = (r*i, ...,fx)- For the given R, 
we calculate the state vector p( a \q)\$> > by performing exactly the spin-isospin algebra 
with the methods developed in refs. [32,44]. The momentum-dependent terms in p^ are 
calculated numerically; for example, 

V i , a *(R) = -^[*(R + 6 i , a )-*(R-6 i , a )] , (33) 

where 5i yCe is a small increment in the ri jCt component of R. The /^-integration is carried 
out with Monte Carlo techniques by sampling a large set of R configurations with the 
Metropolis algorithm. 

The two-body pion and p-meson operators have been constructed from the Argonne 
V14 following the method outlined in ref. [44]. This implies replacing the propagators in 
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eqs. (22a-c) by the Fourier transforms v aT (k) and v tT (k) of the isospin dependent spin-spin 
and tensor components of the interaction model as 



a 2 1 



4m^ /c 2 + m 2 

gjLiv(i + 'v) 2 i 

4m 2 /c 2 + m 2 



K-(fc) = 2^ T (A;) - v (JT (k) (34a) 
V^(fc) = u tT (/c) + v aT (k) (34b) 



The replacements eq.(34) are the ones required for the construction of a two-body elec- 
tromagnetic current operator that satisfies the continuity equation with the interaction 
model [44] . We here apply this replacement to the pair current EM and strangeness charge 
operators as the generalized propagators constructed in this way are then consistent with 
the short-range behavior of the corresponding interaction components. This short-range 
behavior is determined phenomenologically by fitting NN elastic scattering data. An ad- 
ditional justification for using the construction eq.(34) is that it has been shown to lead 
to predictions for the charge and magnetic form factors of the trinucleons [32,42,44], and 
threshold electrodisintegration of the deuteron [45] that are in reasonably good agree- 
ment with the empirical data. The w-meson propagator in the corresponding pair current, 
Eqs. (22b, c), and the p- and 7r-meson propagators in the transition current, Eq.(23), are 
modified by the inclusion of monopole meson-nucleon form factors 

2 A 2 M - m 2 M 
F NNM (k ) = p — ^ 2 M , (35) 

where M is the exchanged meson of mass m M and A M is a cut-off parameter. We use the 
values Ann w = Ann p = Ajvjvtt = 2 GeV, as obtained in boson exchange interaction models 
[46] . It should be emphasized that the contributions due to the vector meson pair currents 
are not significant in the momentum transfer range of interest here. Furthermore, we note 
that in evaluating the contributions due to the vector meson pair currents that include 
the retardation correction, the non-local terms in eq.(22c), namely those proportional to 
p + p', have been neglected. This is justified for the p-meson pair current, since the 
non-local contribution is suppressed by a factor (1 + k p ) 2 (k p = 6.6) with respect to the 
leading term proportional to F^ s \t). This approximation, however, is questionable for the 
o;-meson pair current, since in this case the tensor coupling is small, k u = —0.12. 

IV. Results and Discussion. The results of our VMC calculation are displayed in 

- (s) 

Figs. 2-6. In computing various contributions to F^ and Fqq, we have employed a 
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value of X E = X n = 5.6 to serve as a point of comparison, although is essentially 
a free parameter characterizing the next-to-leading <5 2 -dependence of Ge and is to be 
constrained by experiment. 

Assuming the values p s = —2.12 and p s = —0.2 for the strangeness radius and mag- 
netic moment of the nucleon, we find that the relativistic Darwin-Foldy and spin-orbit 
corrections to the single nucleon operator, and the two-body contributions associated with 
pseudoscalar and vector meson exchanges as well as the pn transition current lead to about 
0.5% decrease (increase in magnitude) of p s [ 4 He], a negligible effect. 

Results for Fq =0 = 2y / 7rF^° and = 2^/ttF^ over a range of momentum-transfer 
are shown in Figs. 2 and 3. Panels 2(a) and 3(a) give the full form factor resulting from the 
one- and two-body currents as well as in the impulse approximation (IA) for comparison. 
Panels 2(b) and 3(b) display individual contributions from the various one- and two-body 
terms. As indicated by the plot in Fig. 2(a) and as noted in previous work [32], the 
inclusion of MEC's significantly improves the degree of agreement with the data on Fq =0 
over a wide range of q as compared with the IA form factor. The difference in behavior 

— (s) 

between F^° and F^ at low-q is dictated by the different values of the corresponding 
nuclear charges: F£=°(0) = AG T E =°(0) = 2 and f£\o) = AG { ° ] (0) = 0. At larger values 
of q, the nuclear EM and strangeness form factors manifest similar structures, having 
their first diffraction mimina and subsequent maxima at essentially the same values of 
momentum transfer. Since the various contributions to F^, =0 (q) are discussed elsewhere 
[32], we focus on F^\q). At low momentum transfer, the nuclear strangeness form factor 

(s) 

is dominated by the single nucleon contribution proportional to G E . In this regime, the 
largest corrections arise from the spin-orbit and p — tv transition currents. At moderate 
values of momentum transfer (q~^2 fm), the largest corrections are due to the pionic pair 
and p — tv transition currents. In arriving at the results shown in this figure, we have 
assumed essentially the Jaffe value for the nucleon's strangeness radius (p s —2) and a 
value of p s = —0.2. Under this assumption of what would be a large magnitude for p s , the 
one-body, pionic, and p — tv transition contributions are of the same order of magnitude 
at the kinematics of the approved CEBAF experiment [24] (q = 3.93 fm -1 ). At this 
point, the p-ir contribution makes up about 20% of the total Fq\ Were we to employ, 
instead, the results of the kaon-loop estimates of the strangeness parameters, p s ps 0.4 and 
p s ~ —0.3 [20], the magnitude of F^ would be an order of magnitude smaller and would 
be dominated by the p-n contribution. 
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We emphasize that the relative importance of the pionic operator is highly dependent 
on one's model for the one-body strangeness form factor, F± s \ which enters the two-body 
operator multiplicatively. In this case, the scale of the two-body operator is set by the 
Dirac one-body strangeness radius, 



Dirac _ ar l 
Ps ~ dT 



Sachs 
s 

T = 



p S s achs + Vs ■ (36) 



In arriving at the results displayed in Fig. 3, we used essentially the pole model value [16] 
for /9^ irac pa —2.4, but not the Q 2 -dependence for F± s \ since the latter is un-realistically 
gentle in light of simple quark counting arguments. Had we used, instead, the results of the 
kaon loop estimate of Ref. [20] , the magnitude of the pionic contribution would have been 
a factor of 20 smaller than the contribution shown in Fig. 3, and the sign would have been 
opposite. Similarly, the one-body IA contribution would be reduced by at least a factor of 
four in magnitude and its sign would also have been opposite than what appears in Fig. 
3. In this case, the p — % transition current would generate the dominant contribution to 
Fq^ at the kinematics of the approved CEBAF experiment [24], while the single nucleon 
strangeness radius would still govern the low-q behavior of the nuclear strangeness form 
factor. 

By way of comparison, we note that the vector meson pair current contribution to 
F^ is negligible at moderate values of q. Although the precise numerical values of their 
contributions depend on one's model for G$ and G^\ as well as on one's choice as to the 
treatment of the retardation term, the overall magnitude of the vector meson pair current 
contribution is sufficiently small so as to render the impact of these model-dependencies 
negligible. 

In Figures 4 and 5, we plot the ratio R s = F^ (q) / F^° (q) = Fq \q) / 'Fq =q '(q) , which 
characterizes the s-quark corrections to the non-strange PV asymmetry (Eqs. (3) and (6)). 
Assuming \G$ /G™\ ~ 1 and \G$ /G 1 [p \ ~ 1> which essentially corresponds to assuming 
the Jaffe values for p s and p s but a more realistic momentum-dependence in the strange 
form factors, we expect a 35% correction to the non-strange asymmetry (the first term on 
the right side of Eq. (6)) at the kinematics of the CEBAF PV 4 He experiment. Fig. 4 
shows the dependence of this correction on the value of p s which, under our form factor 
parameterization (Eqs. (12-13)), sets the scale of contributions generated by G$ '• For 
purposes of illustration, we have assumed magnitudes and relative signs G^/G™ pa —1 
and Gm /G% =0 pa —1 and have taken a positive sign for g pns - The results at low-q imply a 
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negligible dependence of on p s in this regime. For momentum transfers in the vicinity 
of those suggested for the approved CEBAF 4 He experiment, the ratio changes by < 15% 
as Us is varied over a range of values suggested by model calculations [16-21]. Had we 
assumed the kaon loop values for p s and p, s , the overall size of the ratio R s would have 
been nearly ten times smaller, so that in this case the relative impact of any uncertainty 
in p s would be correspondingly enhanced. 

Fig. 5 displays the impact on R s made by the choice of sign of gp-K S , which one cannot 
determine from 0-decay data and symmetry arguments. We illustrate this sensitivity for 
two different models of nucleon strangeness: (A) (p s , p, s ) = (—2.0, —0.2) and (B) (p s , p, s ) = 
(0.0, —0.2). For low momentum-transfer, the impact of this uncertainty in sign is negligible, 
whereas at q ~ 4 fm _1 (the kinematics of Ref. [24]), it corresponds to roughly a ±15% 
uncertainty in the asymmetry. To put the point somewhat differently, even if the strange 
quark vector current matrix elements of the nucleon vanished identically, we would expect 
non-nucleonic strange quarks in the nuclear medium to generate a 15% correction to the 
non-strange P V asymmetry at the kinematics of Ref. [24] . The scale of this effect is well 
below the 40% statistical error projected for the approved CEBAF experiment, assuming 
50% beam polarization (the error is reduced to 28% for 70% beam polarization). Thus, 
for a measurement of A iH ( 4 He) to be sensitive to gpJ(Q 2 ), significantly a longer running 
time and/or higher beam polarization would be required. 

In Fig. 6, we present the significance of a moderate-g A LR ( 4 He) measurement from 
a somewhat different perspective. If one wishes to constrain the various strangeness pa- 

(s) 

rameters p s , p s , and A^ at a level necessary to test model predictions in detail, then a 
combination of experiments using proton and A > 1 targets would be required [1,2,22]. 
As noted in Refs. [1,2,22], a combination of low- and moderate-g PV experiments with 
4 He could potentially constrain (p 8 , \ J! ) more tightly than could a sequence of A LR (ep) 
measurements alone. This conclusion was based on a one-body (IA) calculation and the 
ideal assumption of 100% beam polarization with experimental errors being statistics dom- 
inated. The inclusion of two-body currents does not alter our previous conclusion about 
the possible constraints attainable from a low-q measurement, since the two-body contri- 
bution is negligible in this regime. In Fig. 6, we display the impact of two-body currents on 
constraints attainable at moderate-g. Fig. 6a shows the joint constraints on (p s , p s ) a 10% 
measurement of A iH ( 4 He) could produce, assuming the parameterization of Eqs. (12-13) 
and central values for these parameters given by model (A) discussed above. A similar plot 
for (p s , \$) constraints is given in Fig. 6b, where a central value for (p s , A^) = (—2, A n ) 
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is assumed. The solid and dashed lines give the constraints corresponding to different 
choices as to the sign of g p7rs . We take the difference between these two sets of lines as one 
measure of the theoretical uncertainty associated with our calculation. 

From Fig. 6a, we note that the correlation between p s and p, s is weak. This feature 
follows from the relatively small magnitudes of the vector meson pair current and spin-orbit 
contributions, which carry the strongest dependences on p, s . In the case of Fig. 6b, we 
observe that the moderate-g constraints are modified only slightly from the IA expectation, 

(s) - 

even though many-body currents generate significant contributions to F^q and F^°. The 
reason for the insensitivity of these constraints to the two-body currents can be explained 
in the following manner. First, the pionic corrections are proportional to the Dirac form 
factor 



F [ a) = (1 + r)" 1 [GT + rGTJ , (37) 

where "a" denotes either the isoscalar EM current or strange quark current. At the kine- 
matics of the moderate-g CEBAF PV experiment, one has r ps 0.17 so that rG^ =0 / 'G T E =0 ps 
0.15. In this case, F[ =0 ps G e =0 . Similarly, tGm/Ge ~ Ps/ Ps ~ 0.15, assuming the Jaffe 
values for the strangeness parameters, so that ps G^. Under these assumptions, the 
pionic pair currents give the dominant correction to the IA nuclear form factors, so that 
at r = 0.17 (q 2 = 0.6 (GeV/c) 2 ) one has 

FSo°(q) « (g-s. ||M T ^) [1] +M T= °(?)SLicll S-S-> (38a) 
« GT°(t) <g.s. \\d( q )^+d( q )W\\ g.s.) 



F^(q) ps (g.s. ||M W ( ? )W + M<S%)!Licll g-s-) (38b) 
psGL s) (r)(g.s. \\d( q p+d(q)W\\g.s.) , 

where O(q)^ and 0(q)w are nuclear operators (see, e.g., Eqs. (15) and (A.l)). Hence, 
the ratio R s = F^q/F^q is essentially independent of nuclear matrix elements and is 
governed by the ratio of single nucleon form factors as in the IA case. Thus, the inclusion 
of two-body currents does not seriously change the joint constraints on p s and A^. Some 
changes from the IA results do appear, since neither G^m nor the p-rr constributions are 
completely negligible. In the event that \p s /p s \ >> 0.15, however, this argument would 
break down and our conclusions would have to be modified. 

We also point out that the uncertainty in the sign of the p-ir transition current con- 
tribution does not seriously affect the (p s , X^) constraints, even though the magnitude of 
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this contribution is as large as the experimental uncertainty in A LR assumed in obtaining 
the plots of Fig. 6. To understand why this is the case, consider the following argument. 
If one assumes that all of SA LR translates into an uncertainty 5R S in the extracted value 
of -R s , and if one further assumes that p s and can be varied in a manner consistent 
with this uncertainty 5R S (as we've done in obtaining the lines in Fig. 6b), then one has 

5 4^ = (98) 

A LR 4sin 2 9 W + R S 

or 

SA 

SR S = (4sin 2 9 W + R s ) « 0.1 + 0.1 R s (40) 

A LR 

for SA LR /A LR = 0.1. Since R s changes by only ±0.1 for different choices for the sign of g pns , 
the impact of this choice on the magnitude of 5R S - and, therefore on the joint constraints 
on (p s , Xe) - is an order of magnitude smaller than the impact of the experimental error 
in A LR . 

(s) 

Finally, in Fig. 7 we show the p-meson pair current contribution to i 7 ^ under the 
two different assumptions as to the inclusion of the retardation correction. The curve la- 
belled by was calculated without the retardation correction (Gari-Hyuga convention 
[36]), while the curve labelled by "p" includes it (Riska convention [32,33]). The differ- 
ence between the two should be taken as an estimate of the theoretical uncertainty in 
the treatment of these short-range currents. Fortunately, the scale of the vector meson 
contributions is sufficiently small that the choice of convention has a negligible impact on 
the value of Fq . 

V. Conclusions. We have computed MEC contributions to the 4 He strange quark elas- 
tic form factor, F^ s) = 2^F^(q), using Monte Carlo methods and an accurate varia- 
tional ground state wavefunction. Our results indicate that the nuclear strangeness ra- 
dius, p s [nuc], which governs F@ (q) at low momentum-transfer, is (1) dominated by the 
single nucleon strangeness radius, (2) two orders of magnitude less sensitive to many-body 
strangeness currents, and (3) independent of pionic MEC's - results which essentially con- 
firm our previous conclusions based on the shell model calculation. At moderate values 
of q, such as those corresponding to the approved CEBAF elastic PV 4 He experiment 
[24], we find that F^ generates a 35% correction to the PV asymmetry, assuming that 
\Ge/G e \ & 1 and \G^m /Glp°\ w 1 in this regime and that g p s J(Q 2 ) is correctly given by 
</>-meson dominance. The magnitude of this correction is smaller than the statistical error 
projected for the CEBAF experiment under the most conservative assumptions about beam 
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polarization. In the absence of nucleonic strangeness, non-nucleonic ss pairs would gen- 
erate roughly a 15% correction to the non-strange asymmetry at these kinematics. Thus, 
the scale of the strange-quark contribution to A iH ( 4 He) is still sensitive to the nucleon's 
strangeness electric form factor. In the event that \Ge /G%\ « 1, a more precise 4 He 
PV measurement could probe the p-n strangeness charge operator. Such a measurement 
would be interesting since only a p-ir transition three-current operator has been probed 
in other experiments performed to date [27]. Finally, inclusion of MEC contributions to 
Fq and F^, =0 does not appear to affect noticeably the constraints on the leading and 
next-to-leading <5 2 -dependence of attainable with a medium-g measurement of the 
4 He PV asymmetry. 
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Appendix 

In this Appendix, we provide complete expressions for the two-body charge operators, 
expressions for the low-g forms of the corresponding Coulomb projections, and additional 
details of our shell model calculation of the nuclear strangeness radius. 

Two-body charge operators 

Expressions for the co-ordinate space charge operators can be obtained by Fourier- 
transforming the momentum-space operators in Eq. (22) and (24) and summing over all 
nucleon pairs. The resulting formulae are 



[2] 



piomc 



— % 



9-ITNN 

167rm^ 



(A.l) 



i<j 



X 



'3 



for the pionic current, 



[2] (a) 
V— pair 



9vNN 

167rm^ 



-rayrij 



G^\r)Y,fv{hj)^i - xl)S{Sj - £/) (A.2a) 

i<j 

e *?-Zi (^ . f i:j + a j x q ■ Oi x fij) 



e -m v rij 


(1 +m v nj) 









+ i(l + K v ) 

- e^ Si (q- f i:j + <Ji x q- <jj x f^) 
for the vector meson pair current in the absence of the retardation term and 



[2] (b) 
V— pair 
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(1 + nivTij) 



x [e^' Xj <Tj x q- diX fij - e l $' Xi (Ji x q- dj x fy] 

with the retardation term included. The p-n transition current operator is given by 

(°) t \ 



?>2-Km 2 N m p 



Y^K^i - ~ ■ Tje^TiiJ) (A.3) 



i<j 
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j) = i Fi(fij)Vj ■ (fy x q)&i ■ q- Fi(-ry) Vj • (f^ x q)Bj ■ q 

+ F 2 (fij)Vj ■ (hj x q)hj ■ Oi + F 2 (-r ij )V l ■ (hj x (z)^,- • 
+ F 3 (f tJ )Vj ■ {&i xq) + F 3 (-f lJ )V l ■ {jjj x q) 

+ (1 + Av p )|<Ji • <jj [i {G'i(fy) + Gi(-fij-)} + g • fij {H^fij) - i^i(-ry)}] 

+ q- Biq- ffj [i {G 2 (f lj ) + G 2 {-rij)} + q- hj {H 2 {fij) - H 2 (-fij)}] 
+ hj ■ Wis ■ dj [i {Gsihj) + G 3 {-fij)} + q- fij {H 3 (fij) - H 3 {-r i3 )}\ 
+ hj ■ Btq- ffj [{G 4 (fij) ~ G 5 (-f i: j)} + iq- {H 4 (fij) + H 5 (-f tJ )}] 

+ f^ ■ a 3 q- a l [{G 5 (fy) - G 4 (-fy)} + iq- fij {H 5 (f i3 ) + H 4 (-ri 3 )}]^ 
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F 1 (f ) = -( 9l + 9 f) 



F 2 (r) = h + 



9o 



771 / -*\ 9o 
F 3 {r) = 

r 



Gi(f) 

G2(?) 

Gs(r) = q 2 



1 ( 9o \ 2 
r{ 9l -YJ q 
2tfi L h 



r 2 



i(»-f) + (fc- 



^0 
2 



G 4 (r) = ko — — 
r 

G^ = — + q \92- ^) 



(A.5a) 



(A.5b) 



#i(r1 = -(V^ 
r V r 



H 5 (r) 



r 
9o 

T -92 



r) 



Sh 3g 
ko + — + -^2 



(A.5c) 



23 



and where 



,1/2 

g n (r, q) = / d(3(3 n exp [i(3q- r- Lr] (A.6) 

J -1/2 
,1/2 

h n (f, q) = d(3(3 n L exp [ij3q ■ r - Lr] 

J-l/2 
,1/2 

k n (r, q) = / d(3(3 n L 2 exp [i(3q- f- Lr] 

J -1/2 



with 



We define 



L 2 = hml + m 2 p )+[3(m 2 p -ml) + (l/4-[3 2 )q 2 . (A.7) 



V i = i V i -^ i , (A.8) 

where V i and Vj are gradients acting to the left and right, respectively, on the co-ordinate 
of the i-th nucleon in the wavefunction (and not on the co-ordinates appearing in the 
operators). The isospin tensor T v {i,j) is defined in Eq. (23), the quantities Xi and x[ are 
the co-ordinate of the i-th nucleon in the initial and final state wavefunction, respectively, 
and where the co-ordinates 7% etc. are defined in Eq. (A. 11) below. As elsewhere, the 
superscript "a" denotes either the T = EM current or strangeness current operators. 

Expressions for the pair current operators with hadronic form factors included (finite 
A M ) may be obtained from the above formulae by making the replacement 

(A 2 - m 2 ) d - 

0(m M )^0(m M )-0(A M )+ { M m) 0(A m ) , (A.9) 

zA M ClA M 

where <D{m M ) is any one of the pair current operators in Eqs. (A. 1,2) associated with 
the exchange of a meson having mass m M . Similarly, for the p-n transition current, the 
Coulomb operator in the presence of hadronic form factors arises from making the replace- 
ment 

6{m^,m p ) -> d{m^,m p ) +0(A 7T , A p ) - 6(m^,A p ) - £>(A^,m p ) , (A.10) 



where Oija^.mp) is the operator appearing in Eq. (A. 3). 
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Substituting the above expressions for the charge operators into Eq. (7c) and ex- 
panding the exponentials in powers of q leads to the following expressions for the leading 
(/-dependence of the Coulomb operators: 
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L*Z = -iRij x V fli 



Eij- = ^ + (A. lid) 

— * 

A — » — * 

ij — ® 3 1 

and where 

Z x (i,j) = (m.n/e-" 1 '^ - (m p rij) 2 e- m " r ^ (A.12a) 



Z 2 (i, j) = [3 + 3(m p r lJ ) + (m^) 2 ] e" m ^ (A.12b) 
- [3 + 3(7^7^) + (m^j) 2 ] e~ m ^ 

Z 3 (i,j) =[2 + 2(m p rij) + (m^) 2 + (m p r, 3 f] e ~ m ^ 

-[2 + 2(m n r tJ ) + (m^) 2 + (m^ry) 3 ] e~ m ^ (A. 12c) 

Z 4 (i, j) = [1 + m p r y ] e"™^ - [1 + rn^ry] e~ m ^ . (A.12d) 
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We note in passing that the overall normalization of the p — n operator appearing in 
Eq. (A. 10c) differs by a factor of four from that appearing in Eq. (7) of Ref. [26]. The 
latter, as well as the terms involving gpj in Eqs. (10) and (14) of that work, should be 
multiplied by 1/4. 

Shell Model Calculation 

Use of a simple shell model 4 He ground state allows one to obtain analytic expressions 
for the nuclear strangeneness radius, p s [nuc], which are useful in the interpretation of 
the numerical results obtained with variational ground state wavefunctions. To that end, 
we compute matrix elements of the one- and two-body Coulomb operators in the low-g 
limit. From the expressions in Eqs. (A. 10) and as noted in the main text of the paper, 
the vector meson pair current and p-n transition current Coulomb operators go as q 2 for 
low-g. The two-body pionic operator, in contrast, vanishes as g 4 since ~ g 2 and 
since the operator carries an additional, explicit factor of g 2 . Similarly, the vector meson 
pair current retardation term also enters at 0(q 4 ). Thus, for purposes of computing two- 
body contributions to p s [nuc] , we need only compute matrix elements of the operators in 
Eqs. (A.10b,c). In the limit that the 4 He ground state is a pure S-state, the leading g 2 - 
dependence of the one-body matrix element (Eq. (15)) is given by the one-body form factor 
times the number of nucleons and is independent of details of the nuclear wavefunction. 
The two-body matrix elements, on the other hand, are structure-dependent. An important 
consideration in this respect is the short-range repulsion between nucleons. Since the 
ranges of the p- and w-mesons are commensurate with the radius of the repulsive core 
in the N-N potential, matrix elements of the vector meson exchange operators ought to 
be suppressed. To account for this effect, we compute the two-body shell model matrix 
elements by including a phenomenological correlation function, g(r), in the integral over 
relative co-ordinates: 



/ r z dr u* (r)Ou(r) — > / r 2 dr g(r) u* (r)Ou(r) , (A.13) 
Jo Jo 

where u(r) is the radial wavefunction for the relative motion of two nucleons, r = \x± — x*2\ 
is the relative co-ordinate, and O is an r-dependent two-body operator. Following the 
approach of Ref. [34] , we take the correlation function to have the form 

g (r) = C \l - e~ r2/d2 ] , (A.14) 
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where the constant C is determined by the requirement that the wavefunction be normal- 
ized. A fit to the nuclear matter correlation function of Ref. [47] gives d = 0.84 fm. With 
this form for g(r), the un-correlated two-body matrix elements are modified as 

M(6) -> C [M(b) - M(6eff)] , (A.15) 

where M(6) is the un-correlated two-body matrix element computed using an oscillator 
parameter 6, where the effective oscillator parameter is given by 



and where 



freff 
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1 + 2i l 



-1/2 
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(A.16) 



(A.17) 



of Fqq the expression given in Eq. (26). The nuclear Ai^ are given by 



Ai = 2/v^ , 



A 



2 a 



V=p,u> \ 



24tt 2 I \m N J (m v b) 
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(A.18b) 
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(m p 6) 2 - (m^ft) 1 



(A.18c) 



x [(m„b) 2 I(m„b) - (m p b) 2 I(m p b)] , 



where 



I(mb) = 1- \/ ^(mb) exp[i(m6) 2 ]erfc (^=^J , 



(A.19) 



and where M v ,2 are spin-isospin matrix elements and g VNN is the vector meson-nucleon 
coupling. For b = 1.2 fm, one has m p 6 >> 1, so that the function in Eq. (A. 18c) may be 
expanded as 



(m p b) 2 I(m p b) = 1 - 



+ 



(A.20) 



(m p bY 

A similar expansion in powers of l/(m v b) has been used in arriving at the expression in 
Eq. (A. 18b), where the + ■ • • indicates contributions from terms higher order in l/(m v b). 
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Figure Captions 



Fig. 1. Two-nucleon (N and N') meson exchange current (MEC) contributions to nuclear 
matrix elements of the isoscalar EM and strange-quark vector currents. "Pair current" 
processes, shown in (a,b), arise from the negative energy pole in the intermediate-state, 
single nucleon propagator. "Transition current" contributions (c) are generated by matrix 
element of the current operators (indicated by the ®) between mesonic states (M' and 
M). As explained in the text, mesonic matrix elements of J^ M (T = 0) and s7 M s vanish 
when M' = M. 

Fig. 2. 4 He elastic charge form factor, F(j =0 (q). Panel (a) gives the absolute value of the 
form factor. Circles indicate experimental values. Dashed curve gives theoretical prediction 
in the impulse approximation (IA) while the solid curve results from the inclusion of two- 
body currents (IA+MEC). Panel (b) shows individual one- and two-body contributions 
to F^ =0 (q). One-body contribution is indicatd by solid curve (IA). Dashed curves give 
contributions from the pionic (circles), p- meson (squares), and w-meson (asterisks) pair 
currents as well as the p-n "transition current (triangles). Short dashed curve indicates 
the spin-orbit contribution. Only the absolute value of each contribution is plotted, and the 
signs are indicated in parenthesis. Vector meson pair current contributions are computed 
including the retardation correction to the charge operator (Eqs. 22c and A. 2b). 

Fig. 3. Same as Fig. 2, but for the elastic strangeness charge form factor of 4 He, 
Fg(q). In this case, only theoretical predictions are shown, since no measurements have 
as yet been made. Computations were carried out using (p s ,p s ) = (—2.0,-0.2), which 
correspond roughly to the pole model predictions for these parameters [16], and a Galster- 
like parameterization for the (/-dependence of the one-body strangeness form factors. A 
positive sign for g piTS was also assumed. 

Fig. 4. Elastic strangeness to EM charge form factor ratio, F^\q) / Fq =0 (q) for different 
values of nucleon strangeness magnetic moment, p s , and fixed strangeness radius, p s . 

Fig. 5. Same as Fig. 4, but for fixed p s and variable p s . In each case, results using the 
impulse approximation for F^ are shown along with results including two-body currents 
for two different choices of sign on (7 p7rs . Panel (a) assumes a large negative value for p s , 
while panel (b) gives the ratio for vanishing nucleon strangeness radius. 

Fig. 6. Prospective constraints on single nucleon strangeness parameters from a 10% 
measurement of the 4 He elastic, PV asymmetry at Q 2 = —0.6 (GeV/c) 2 . Solid (dashed) 
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lines give the band of allowed values for positive (negative) sign on g p7rs . In panel (a), 

is assumed fixed, while in panel (b), /j, s is held constant. In both panels, central values of 

(p s , fi s , A^) = (—2.0, —0.2, 5.6) are assumed for purposes of illustration. 

Fig. 7. p-meson pair current contribution to (q) computed including the retardation 
correction and omitting it ("#" subscript). Only absolute value is plotted, while sign is 
indicated in parenthesis. 
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